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CONVERGENCE OF A SEMI-LAGRANGIAN SCHEME FOR THE BGK 
MODEL OF THE BOLTZMANN EQUATION 

GIOVANNI RUSSO, PIETRO SANTAGATI, AND SEOK-BAE YUN 

Abstract. Recently, a new class of semi-Lagrangian methods for the BGK model of the 
^vj . Boltzmann equation has been introduced [5] 1171 [T5] . These methods work in a satisfactory 

way either in rarefied or fluid regime. Moreover, because of the semi-Lagrangian feature, 
r^ ' the stability property is not restricted by the CFL condition. These aspects make them 

very attractive for practical applications. In this paper, we investigate the convergence 
vQ ' properties of the method and prove that the discrete solution of the scheme converges in 

a weighted L^ norm to the unique smooth solution by deriving an explicit error estimate. 
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1. Introduction 

In the kinetic theory of gases, the dynamics of a non-ionized monatomic rarefied gas 
system is described by the celebrated Boltzmann equation. But numerical approximation of 
the Boltzmann dynamics is a formidable challenge due mainly to the complicated structure 
of the collision operator. Many good numerical techniques have been developed to this end, 
but often they lead to time consuming computations. 

To circumvent these difficulties, Bhatnagar, Gross and Krook [1], and independently 
Welander |20] . proposed a simplified model for the Boltzmann equation where the collision 
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^^^^ % + ^-^-f= ^(-^(Z) - /)' (^' ^'*) 6 T'^ X M'^ 



operator was replaced by a relaxation operator: 
+ vV,J = -{M{f\ 

K 

f{x,v,0) = fo{x,v). 

Here T*^ denotes the d-dimensional torus and — is the collision frequency. Although the 

k 
collision frequency takes various forms depending on hypotheses imposed in the derivation 

of the model [3l [HI [191 (HI (HI [23] , we assume in this paper that it is a fixed constant for 

simplicity. Ai denotes the local Maxwellian constructed from the velocity moments of the 

distribution function /: 

p{x,t) ( |u-C/(x,t)p' 



M(f){x,v,t) = , ^^ ' ' ^exp -J ^; '/' 

^^'^ ' ' . (iTTTixf.)^ V 2Tx.t) J 



^{2TrTix,t))d "V 2T{x,t) 



where 



Pix,t) = / f{x,v,t)dv, 
p{x,t)U{x,t) = / f{x,v,t)vdv, 
dp{x,t)T{x,t) = / f{x,v,t)\v-U{x,t)\^dv. 



The BGK model (jl.ip is computationally less expensive than the Boltzmann equation 
since it is sufficient to update the macroscopic fields in each time step. On the other hand, it 
provides qualitatively correct solutions for the macroscopic moments in fluid regime. These 
two aspects, namely, the relatively low computational cost and the correct description of 
hydrodynamic limit, explain the interest in the BGK model and its variations over the last 
decades. It also shares important features with the original Boltzmann equation, such as 
the conservation laws and the dissipation of entropy: 

1 \ . / I 



(1.2) / M{f) I V \ = I f\ V \dv 

and 



\v\ / \ \v\ 



(1.3) jiMif)-f)\ogfdv<0. 

(II. 2p leads to a sy 
■ f fdv + V,- fvfdv = 0, 



The local conservation laws (II. 2p leads to a system of hydrodynamic type equations. 

d_ 
dt 

(1.4) — I fvdv + V^- I v(^vfdv = 0, 

d 

lit 

There are extensive literatures on various topics of the BGK model. For mathematical 
analysis, we refer to [51 [9l [121 [IH [HI [HI [12]- For numerical computations, since there 
are too many of them, we do not attempt to present a complete set of references. See 
[6l [m [131 [3 [S] and references therein. 



■ I f\v\^dv + V^- f v\v\'^fdv = 0. 
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Recently, a semi-Lagrangian scheme was proposed and tested successfully for various flow 
problems arising in gas dynamics [H \T7\ I18| . The first order version of the method can be 
written down as follows: 



where M^j denotes the local Maxwellian defined by 



~ n" / \v- — [/" P 



«J 



The precise definitions of each terms will be given in later sections. The main feature of the 
scheme is that even though the relaxation operator is treated implicitly, computations can 
be performed explicitly by exploiting the approximate conservation laws in a very clever 
way. (See section 2). Therefore, (jl.Sp enjoys the stability property of implicit schemes 
and the low computational cost of explicit schemes at the same time. Moreover, the semi- 
Lagrangian treatment of the transport part enables one to perform the computation over a 
wide range of CFL numbers. In this paper, we study the convergence issue of this scheme 
and derive an explicit estimate of the convergence rate measured in a weighted L^ space. 
As far as we know, this seems to be the first result on the strong convergence of a fully 
discretized scheme for nonlinear collisional kinetic equations. 

This paper, after introduction is organized as follows. In section 2, we describe the 
numerical method considered in this paper. In section 3, we recall relevant existence results. 
In section 4, we present our main result. Then several essential estimates to be used in later 
sections are presented in section 5. In section 6, we derive a consistent form and obtain 
error estimates of the remainder terms. Finally, in section 7, we combine these elements to 
prove the main theorem. 

2. Description of the numerical scheme 

For simplicity, we consider one dimensional problem in space and velocity. We assume 
constant time step Ai with final time Tj and uniform grid in space and velocity with mesh 
spacing Ax, A-y respectively. We denote the grid points as follows: 

t" = nAi, n = l,...,7Vt, 

(2.6) Xi = i/\x, i = l,...,Nx, {modi), 

Vj=jAv, J = -iV^,,..,0, ..,iV^,, 

where NfAt = Tf, N^Ax = 1 and N^Av = R. We also denote the approximate solution 
of f{xi,Vj,t) by fl^j j^. To describe the numerical scheme more succinctly, we introduce the 
following convenient notation. First, we define x(i,j) = Xi — Atvj. We also set s = s{i,j) 
to be the index of the spatial node such that x{i,j) G [xs,Xs+i)- 

Definition 2.1. We define the reconstructed distribution function ff'j ^ as 

fr. r,\ 7n _ ^V-^J) ~ ^s j-n , Xs+l — X[l,J) 

\^-') Ji,j,R - -^ Js+l,j,R H -^ Is,j,R- 

Note that /"• ^ is the linear reconstruction of /(x — f At, -y, nAt). The numerical scheme 
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■fi^jM 


f/j ; ; ; pi 
Jsj\ : 1 Js-\-ij 




■^.5 






Xs+l 





o- 
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Figure 1. Characteristics diagram for positive velocity grid node. 



we consider in this paper can now be stated as follows. 
(2.8) 



"" f^R+T^.^UfR)^ 



Ji,j,R K + At '^'J'" ' K + At 



where 



fi,j,R= {f0'^{\v\<R}){xi,Vj) 

— fl 



■n \N 



C^^T^r) 



\^J ^i,R\ 
'^-'-i,R 



and 



(2.9) 



K 



i 
j 



We now explain briefly how (j2.8p is derived. The numerical scheme for (jl.ip is based on 
the following characteristic formulation of the problem: 



(2.10) 



f = -(-M(/)-/), 
at K 

dx 

—-= V. 
dt 
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Then the time evolution of fij{t) = f{xi, Vj,t) along the characteristic line in the time step 
[tn,tn+i] is presented as 

(2.11) f '' 

dxi 

—— = Vj. 

dt ^ 
We then discretize ()2.1ip implicitly to obtain 

(2.12) "^' ^^ '' = -(-^-ri/r') - /^')- 

At this step, (I2.12P seems to be time consuming, and this is where the clever trick kicks 
in. We first observe that conservation of mass, momentum and energy gives for (f){vj) = 

l,Vj, 2\'Vj\ 



0, 



which implies 



E/^'<^(^^)^E/m-'^(^^- 



This in turn gives 

(2.13) a^i:;h/i7')^-^m(S)- 

We then substitute the above approximation (I2.13P into (j2.12p to obtain 

(2.14) -^"^^r^ = -^^um - fit')- 

Note that the implicit scheme (j2.12p now can be calculated explicitly. We then collect 
relevant terms together to obtain (|2.8p . For more details about ()2.8p , we refer to \n\ [T8] . 

2.1. Extension of the scheme. In this section, we extend the discrete distribution func- 
tion {fij}i,j to the whole numerical domain T^; x M^ and reformulate (j2.8p in accordance 
with the extension. This allows us to treat the discrete numerical solution and the exact 
solution in the same framework. First we introduce 



i 
C2{v) = ^Vi X^^^^^<^^^ y 



Here Xa denotes the usual characteristic function and we used the following convenient 
notation: 

^i±i ^''^^^^ («±2^ 
We now define 



Axiix.v) = Xa^^{x,v), 
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where 

^Xi = {ix,v)\ Cl(x) = Xi}, 

^5 = {ix,v)\C2{v)=vj}n{{x,v)\\v\<R}. 

For given sequences {oj} and {hj} defined on grid nodes, we define the following extension 
operators: 

Ex{ai){x,v) = ^ (^ ' flj+i + -!^ aijAx,{x,v), 

id 

Now the approximate distribution function can be extended to the whole numerical domain 
as follows. 



f^{x,v) = Eifr^^^^)ix,v) 

Ax ■''+''^ ' Ax 



(2-15) = E(^/m. + ^^^^/M>5(->-)- 



«j 



Note that /|j(x, u ) is piecewise constant in the velocity domain and piecewise linear in the 
spatial domain. Using this, we define the macroscopic fields and the local Maxwellian as 
follows: 



PRix,t) = / f^ix,v,t)dv, 

(2.16) plix,t)U^ix,t) = j f^ix,v,t)C2iv)dv, 

dplix,t)T^ix,t) = J f^ix,v,t)\C2{v)-U^ix,t)\^dv 
and 

(2.17) A,»,/S),„).^JftLexp^ lC.i^)-UM^)? 



v'(2Trs(x) 'V 2,rs(i) 

We also define the reconstructed distribution function f^ as 

Then the corresponding macroscopic fields and the local Maxwellian are defined analogously: 



Pjtix,t) = / fR{x,v,t)dv, 
(2.18) ^(x,t)f7^(x,t) = j f^ix,v,t)C2iv)dv, 

drRix,t)f^{x,t) = [f^ix,v,t)\C2iv)-U^{x,t)\'dv 
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and 



(2.19) ■M"(/g)(x,t;)= /^^^^ :exp( 

/ c\ mil / \ ^ 



pUx) _ ( \C2{v)-Ul{x)^' 



(27rT^(x) ^ 27rr^(x) 

2.2. Consistency. The following series of lemmas show that the preceding definitions are 
reasonable. 

Lemma 2.1. /]j, f^ are periodic functions with period 1: 

(2.20) 

Proof. This follows directly from the periodicity of the spatial domain. D 

Lemma 2.2. For each Xi and Vj, the following consistency properties hold. 

flii^ii'Vj) = fi',j,Rj 

(2.21) 

fRy^i^Vj) = fi',j,R- 

Proof. We recall the definition (|2.15p of f]i{xi,Vj) to see 

fR\Xi,Vj) = /_^y ^- fm+l,£~^ ^Z fm,ij-^m,e{^i^''^j) 
m,l. 
^^ ~ Xj ^ ^t+1 ~ Xj „„ 

rn 

The second statement follows in a similar way. D 

Lemma 2.3. For macroscopic fields, we have 

PR{Xi,t) = Pj^j^, 

(2.22) pl{x„tmix„t) = plj,Ur,n, 

PR{Xi,t)Tj:^{Xi,t) = Pi^^T^pf 

and 

p]i{xi,t) = p^ji, 

(2.23) rRix^,tmx^,t)=p^^JiU^,R, 

p]{{xi,t)T^{xi,t) = Pi^nT^R. 
Proof. We prove the third identity. We observe 
dpl{x,)mx,,t) 

f^{x„v)\C2{\v\)-U^{xi,t)\''dv 



m,£ 

/E (^^^"+M + ^^^^fr,i)\C2{\v\) - U^{x.,trAg,ix.,v)dv 
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Y,fr,i\C2i\v\) - U^{xut)\' A^,ix,,v)dv 

e. 

Other identities can be proved in a similar manner. D 

Lemma 2.4. The following consistency properties hold for local Maxwellians. 

M-if^){x,,v,)=Ml^ifr,,R), 
(2.24) 

M^{f^){x,,v,)=Ml^{fr,,R). 

Proof. This follows directly from the definition (I2.19P and Lemma 12.31 D 

Our main scheme (12.81) can now be restated as follows 



m+l ^ ~^' - ^* 



(2-25) f^2k = ^TA^/«(^-^^) + ^TA^-^"(/^)(^-^^)- 

Applying the extension operator to (j2.25p once more, we obtain the following reformulation 
of (Ei 



Theorem 2.1. The discrete scheme \2. ^ can he recast in the following form: 

f'^+\x,v) = -^^Jl{x,v) + -^^E{M'^(h))[x,v), 
(2.26) JR ^ ^ > K + At^^ ' K + At ^ \JR'^y ^' 

fR{x,v) = E{foX{\^\<R}){x,v). 

where we used a slightly abbreviated notation for brevity: 

EiM''{f^)){x,v) = EiM''{f^)ix„Vj)){x,v). 

2.3. Notation. Before we proceed to the next section, we set some notational conventions. 

• C denotes generic constants. 

• Cx^y,.. denotes generic constants that depend on x, y,... but not exclusively. 

• We use the following convention for the L^ norm with polynomial weight and the 
L°° norm. 

ii/wiili = / fix,v,t)ii + \v\rdxdv, 

||/(t)||^oo = sup\f{x,v,t)\, 

x,v 

||/(t)|U^ = sup|/(x,t)|. 

X 

• We introduce the following notation for weighted -L°°-Sobolev norms for smooth or 
approximate solutions. 

N^ifm = sup\ii + \v\rf{x,v)\, 

x,v 

<(/)W = J]sup|(l + b|)«9,"/(x,^)|, 

\a\=k 

Xim) = E ^^p\{i + \v\ro^d^,f{x,v)\, 
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and 



N^imit) = sup\{l + \ve\yfUxm,Vi: 



ml 



N if^m = sup {l + \v,\) 



m,£ 



, fR{Xm+l,Vi) - f^{Xm,Vi) 



Ax 



For simplicity, we set 



N,{fm = iVO(/)(t) + iVj(/)(t), 
N,{f){t) = iVO(/)(t) + iV;(/)(t). 

Remark 2.1. Note that we have deliberately distinguished Nq from Nq. This simplifies 
many computations in later sections. 



3. Existence and uniqueness of smooth solutions 

In this section, we recall relevant existence results of (jl.ip . The existence and uniqueness 
was first obtained in [15] and the regularity was investigated in [9]. The following theorem 
is a slight simplification of the corresponding results in [9l [15], which is enough for our 
purpose. For the proof, we refer to [9[[T5]. 

Theorem 3.1. [3 [15] Suppose /o > 0, /o € L^{T x M). Suppose further that there exists 
constants C\ and C2 such that 

I fo(x- vt, v)dv > Ci > 0, 
(3.27) V 

Nqifo) < C2. 

for p > 5. Then there exists a positive number Tf and a unique solution of the BGK model 
\1.1\) such that 

(1) Nq-norm of f is uniformly bounded in [0, Tj].- 

Nq{f){t)<CT, forte[0,Tf], 

(2) Macroscopic fields satisfy the following estimates: 

Or fO\„CnTf 



Lr + mmLT + \\Tit)\\L^ < CqN:^if'n)e'^"'f 



-C,Tf 



(3.28) p{x,t)>Cqe 

T{t) > Cqe-^^^f > 0. 

In what follows, we list some of the important estimates satisfied by the smooth solutions. 
Readers are referred to |15] for the proof. 
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Lemma 3.1. [15J Under the assumptions of the previous theorem, the following estimates 
hold. 

p{T+\U\^)'^<C,Ng{f) {q>d + 2), 

(^•2^) ^-—I^<C,N,{f) {q<d), 

(T« + |f/]J|2) — 

^' ' <C,N,{f) {q>l). 



4. Main result 

We are now in a place to state our main result. 

Theorem 4.1. Let f be a smooth solution corresponding to a nonnegative initial datum 
/o satisfying the hypotheses of theorem I3.il Let f^ be the approximate solution constructed 
iteratively by ()2.26p . Suppose that the time step is bounded in the sense that: 

(4.30) At < max{-, k} 

and the size of spatial and velocity meshes satisfies the following smallness assumption: 

(4.31) Ax + Av<lM^^^^^^Il^^^, 

2N,{f,){2 + Tf) 

where foR = foX\v\<R- Then we have 



Remark 4.1. L Note that ^4-30 ) is not a smallness condition. Therefore, Theorem \4.1 



shows that our scheme works well even for large time steps. 

2. The condition ^.30 ), ^.31 ) on the mesh size is introduced to derive the lower bound 



estimates on discrete macroscopic fields. See lemma 15.51 

3. For convergence, we need to set the size of Ax and Av to be comparable with At. For 

example, if we set Ax = Av = (Ai)^"^"*, then we have 

ll/(-,-,r/)-/^1lL^<C(r,g,/o)((At)- + At + ^^^-l^). 

4- For high-order methods |17l I18j . we expect to obtain a high-order error estimate of the 
following form: 

\\f{;;Tj)- /^- II,. < c{iAxY + (A.)- + ^^^^' +f ^^"^ + (At)- + ^^-^ ) , 

which we leave for the future research. 
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5. Basic estimates 

In this section, we present several estimates which will be crucial in later sections. Unless 
it is necessary, we do not restrict our argument to one dimensional problems and present 
the result in general d-dimension. 

Lemma 5.1. The reconstruction procedure does not increase the Nq-norm of the discrete 
distribution function: 

NM) < Nqifl). 
Proof, (i) The estimate of N^(f^): We observe from the definition of f^ 



N'qUR) 



SUp\f^{Xi,Vj){l + \Vj\)'^\ 



sup 



Z^^. ■'s(i,j)+l,j,R 



+ ^ Js{i,JUR [^ + ml)' 



Ax 



< 



< 



sup I max{/;(^_^.)^^. ^, /r(ij)+i,,- ^}(1 + \vj\y 



«j 



sup|/;^.^(l + b,|)'i 



«J 



= Nqifi^)- 
(ii) The estimate of Ng{f^): We first define 



2^(^ + l)i) ~ Xs{i+l,j),j _ x{i,j) — Xs[i^j)^j 



Ax 



Ax 



which gives 



l-a 



,(i+ij)+ij - x{i + 1, j) _ Xs(jj)+ij - x(i, j) 



Ax 



Ax 



Therefore, we have from the definition of /]j 

jR{Xi+l,Vj)-jR{Xi,V.j) 



Ax 



^fs{i+l,j)+l,jR + (1 — o)fs{s+l,j),jR~\ \0'fs{i,j)+l,jR + (1 — 0')f! 



Ax 

'fs{i+l,j)+l,jR - fs{i,j)+l,jR 



's(s,j),jR^ 



<a 
< 



Ax 

fs{i+l,j)+l,jR - fs{i,j)+l,jR 



< 



Ax 
aiV,i(/g) (l-aK(/g) 
(1 + 1^,1)'? (1 + l^il)'? 



+ (1 - a) 
+ (1 - a) 



Ax 

fs{s+l,j),jR - fs{s,j),jR 

Ax 

fs(s+l,j),jR - fsis,j),jR 



Ax 



Hence we have 



N^if^t) 



= sup 

< N]{f^). 



(I _|_ \y |\'? '^fi(^'+l'^j) fR\Xi,Vj) 



Ax 
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We then combine the above two estimates to obtain 

NM) = N,{Jl) + Nl{Jl) 

u 

Lemma 5.2. Suppose Nq{f^) < oo with q > d+2. Then the following estimates hold. 

Pi 



q — d 



- < CgNoif^), 



PRin + \U^\')^<C,N,if^) {q>d + 2), 
(5.32) pl 



d—q 

R^\'^R\ '^~ 



PrPrI'^^'^ 



<C,N,{f^) {q<d), 



[{t^ + \u-\^)t:^]-2 



T<CgNgif^) iq>l), 



and 



rR{TR+m')'^<C,N,{f^) iq>d + 2), 
(5.33) fPn 



^ ^ -^<C,N,{f^) iq<d), 



[(r« + |c/]j|2)r]j] 

Proof. We only prove (j5.33p . (j5.32p can be proved in a similar manner, 
(i) The estimate of -^: We have from (IZTHI) 



Pr 



j fl{x,v)dv 

k I _ Jr\C2{v) - mx)\'dv + / _ f^dv 



1 



We take D = ( cn'(A) ) ^ *° obtain 



~n 

^ < C,N,{fl) 



(TR)^ 
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This proves (fO^ . 

(ii) The estimate oi p%{f^ + lU'^l"^)^ : We note that 

fP^{df^ + \U^\^) = jJl\C2{v)-Ul{xtdv + pl{x)Ul{x) 

rR\C2{v)?dv 

C2{v)\>D \'~'2{V)\'J J\C2(v)\<D 



J\C2(v)\>D \'~^2[V)\'' J\C2(v)\<R 



'\C2iv)\>D \'-^2\u)\-^ J\C2iv)\<R 

< C.^^^ + rnD' 



We take 

D 



nMj\—^ 



f^R 



to complete the proof. 

(iii) The estimate of —^ z^ — d^ ■ Note that 



^1 



'fl^l'^H 



PR = [ fldv 

= I Jldv + f Jldv 

J\C2{v)\<D J\C2{v)\>D 

J\C2{v)\<D \<-^2[v)\'i U^ J\C2{v)\>D 

We take D'^'i^^ = p{df^ + |f7^| 2)^-^4+2 to see 

p% < c,N,if^)^^^[p{dn + m\')]'-^^'- 

(iv) The estimate of — r-^ — 5 j- : For o > 1, we have by Holder inequality, 

[(T]J+|l/SP)]^ 



rRm\ < j mc2{v)\dv 



< / . f^\C2{v)\dv+ __ fR\C2{v)\dv 

J \C2{v)~U'^\<D J\C2{v)-U'^\>D 

< {PrY^UI' _ n\C2{vWdv)' + ^ [ \C2{v)-U^\\C2{vmdv 

< C{rR)'-'N,{f^)'^D'^+l(^l\C,(v)\'f^dvf[l\C2{v)-U^^^^ 

< c{rnf~-.N,{rn)hD-. + ^{df^ + m?)'^fi\ 
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We maximize the estimate by taking 

to obtain the desired result. D 

The next lemma shows that the Nq- norm of the discrete local Maxwellian can be con- 
trolled by the Ng- norm of the approximate distribution function. 

Lemma 5.3. Suppose Nq{f'^) < oo with q> d+2, then we have 

NqiM^UD) < CqNqifl). 

Proof. (I) The estimate of iV°(7W"(/g)) : We observe that 

(1 + \C2{vWM-{rR) < Cq{l + \U^\^ + IC2W - U^HM'^if^) 
(5.34) <c (-^ + o"^^^ + oVT^)^ 

By (I5:32]l . we have 

(Q-d.) 



(5.35) -^, plin)^ < CqNqifl). 



The estimate of Op — ^-3- is more involved. We divide it into the following two cases. 



Case 1: |C/^| > (T^)2. 
We have from ([02|) 



\TTn\q \TTn\d+q 

^n \^r\ < n \^r\ 

Pr , rr ^ Pr. tt: :t 



{T^)-2 p^nr^. 



(5.36) Pim'^^'i 



<a 



\Tr^+m¥iTRf^ 



1 

' R\ -y R^ 

We apply ()5.32p to see 



Case 2: \U^\ < (T^')2 



< CqNqif^). 



\TTiT-\1 



(5.37) 






'R}' 



< CqNqif^). 

We combine the Case I and Case II to obtain 

(5.38) pl^ < CqN'qUl). 

We then substitute ()5.35p . ()5.38p into (j5.34p to complete the proof. 
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(II) The estimate of N^{A4{f^)) : We have by Taylor's theorem 
(5.39) 



n 
Pi+1,R 



{2^Tl\ 



:exp 



/_ r^ 



1^-' '-^i+l,R\ 



i+l,R) 



2T 



i+l,R 



Pi,R 



= {f^+i,R - pIr)^-^^ + {u?+i,R - ur,R) 



(2vr7;'^W 

da 



: exp 



U,. _ r/" |2 

ryrpn 

'^■''i,R 



where 



dT ' 



dX 
for some < ^ < 1. We recall that 



dX 



e/iVi,,+(i-e)/", 



dp 



'n\d 



: exp 



i2nT^) 



\v, - U\ 



n\2 

r\ 



2T^ 



to get 






1 



< 



i27rTJ^r 
1 



: exp 



: exp 



2T^ 



Tjn\2 

^r\ 



2T 



(1 + hir 



'^ _i_ U,. _ r/"!"? -I- IT"/""! 



Other estimates can be obtained similarly as follows: 



dU 



{i + \vj\r < c,,T, 



i? 



^(1 + hl)^ < C7, 



g,T- 



R 



(5.40) 



We substitute the above estimates into (j5.39p to see 

< C,{\p^+^,R - pIr\ + \Ur+,,R - Uni,R\ + |7;"+i,R - Tj^j, 
We now estimate each terms separately. First we observe 

\Pi+l,R ~ Pi,R\ = Z^ifi+l,j,R ~ fi,j,R)^'' 



(5.41) 



sup 



-(i + KDIEti 



AxAf 



fn fn 

Ji+1J,R -liJ^R, 

„-, — a;; — ^-I'^^-inz^d + Ki). 



< CX(/]^)Ax. 
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Similarly, we have 

(5.42) \ur+i,R-u::n\ < c.niudax, 

(5.43) \Tt,i^R-T^,R\ < C.NlifD/Xx. 

We substitute (lOTT) . (I5:i2]l . (lOSjl into (fSliOD to obtain 

< c-.iviajj). 

We now establish the stability estimate for the scheme (j2.26p . 
Lemma 5.4. Suppose Nq{f^ < oo with q> d+2. Then we have 



D 



Proof. We take Ng norm on both sides of (j2.26p and apply Lemma 15.11 and Lemma 15.31 to 

At 



'9 

see 



NM) < ;7^^.(/r^)+^^^.(^(/r)) 

K + At ^^■'^ ' 

( g.-l)At x , 



< (1+ ,, ;. )^.(/: 

Iterating the above inequality, we obtain 



-.</«) ^ (>-^^l^)"w«) 

(Cq-l)iVAt 

where we used (1 + x)" < e"^ and nAt < NtAt = Tj. U 

Lemma 5.5. Let q> d + 2 and Ngi^f^) < oo. Suppose that the time step is bounded in the 
sense that 

(5.44) At < max{-,K} 

and the mesh size for spacial and velocity nodes satisfies the following smallness condition: 

(5.45) Ax + A.<LM^^^^Z>1^. 

2N,{f^){2 + Tf) 

Then the following estimates holds for approximate macroscopic fields. 

Wpimk^ + \\Um\\L¥ + WT^Wl^ < CgNM)e('^Tf, 

(5.46) p^(x,t) > C7ge-^«^/, 
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C,Tf 



T^{x,t) >Cqe-^''^f >0 
Proof. Note that we have from Lemma 15.41 

pm = lfR{x,v)dv 






< C,e'''>''iN,U\ 



R) 



To proceed to the estimates for V^ and T]j, we need to estabhsh the lower bound for p^{x, t) 
and T'^{x,t) first. Note that we have from (j2.26p 

PRix)= I mx,v)dv 



(5.47) >___//-i(^,„)d^ 



K 
K~At 



K 



In the last line, we used 
f'^-\x,v)dv 



Ej / f'^-'ixi-C2{v)At,v)dv 



= E { (^ ( E / A^i.-^^./-) + ^^ ( E / ^.'^^^. (-' -H 
= E{(^(Ea>i,H + ^^^^(E^HK(^^ 

= ^^(E&«H 

3 
T? f\^ f ^i-Vj^t-Xs ^n-l , Xs+1 - {Xj - VjAt) ^^^^ -^ 

= ^4Z^I A^ ^^+1'^-.^ + A^ ^^'^ )^V 

3 
P f f\^f Xi-C2{v)^t-Xs ^n-l , Xs+1 - {Xj - C2iv)At) ^^_i\ ARr„.\M, 

= ^4yZ^l A^ /^+i,.,ii + A^ 4. JA,W^^ 

= E.,(^j rR-\xi - C2{v)At,v)dv). 
On the other hand, we note from the definition of E^ that 

^- ( / fR~\^i - C2{v)At, v)dv) > inf f fl{x, - C2{v)At, v, t)dv. 
This gives from (|5.47p 



pI{x) > -—j-inf / f^{x, - C2{v)At,v,t)dv. 
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We iterate the above lower bound estimate to obtain 



Pr{^) > 



(5.48) 



K + At 

K \ " 



inf / f%{xi - C2{v)nAt,v)dt 
K + AtJ '''fjfR^^-C2{v)Tf,v)dv. 



Then we employ the following elementary inequality 

(1 - xY > e"^""^ (0 < X < 1 - e) for < e < 1, 
to see 



K \'n. 



K + At 



1- 



At \n 



K + At 

SnAt 



> e «+^* 



8T 



/ 



> e '^+^* . 



Here we used the fact that the boundedness asuumption on At implies 

1 

2' 



At • fl T 

< 1 — mm |— , KJ. 



Hence (I5.48P gives 
(5.49) 



K + At 



pI{x) > e-^+^^ inf / /)^(x, - C2{v)Tf,v)dv. 



Now we should replace the estimate of f^ with the estimate of /q/j. We first observe that 
the difference between /^ and /or can be estimated as follows: 



fm{xi -VjTf,Vj) 

Xj Vjlj Xg 



< 



Ax 

X i (J 4 J. f X i 



f0R{Xs+l,Vj) + 



Xs+1 - {xi - VjTf) 



Ax ■'^+i'^'^ 



J s 



f^{xi-VjTf,Vj) + 



+ 

dfo 
dx 



^' ' Ax 

Xs+l - {Xi - VjTf) g 



f0R{Xs,Vj) + —-{xg,Vj) 



[Xg,Vj 



<f^{xi-v,Tj,v,) + Ngifo)- 



Ax 
Ax 
Ax 



fs,j,R+ 5^(^^'^i) 



dx 

Ax 



Ax 



{l + lv.il)"' 



where xq £ [xg, Xs+i) and s = s{i,j) denotes the index of spatial mesh grid such that 
Xi — VjTf £ [xs, Xs+i)- We then sum over i to obtain 



^ foR{xi - VjTf, Vj)Av 



(5.50) 



<E/^(^^-^i^/'^j)^^ + E^5(/o)7Y 



AxAf 



(1 + b,- 



/ flixi - C{v)Tf,v)dv + N,{fo)A3 
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On the other hand, apphcation of Taylor's theorem on each interval [v-i,v-,i_) gives 



2 -J ' 2 



/ fonixi - vTf, v)dv - ^ fm{xi - VjTf, Vj)Ar 



(5.51) ^ _i 

<{l + Tf)iAx + Av)NlifoR) 

<{l + Tf){Ax + Av)Ng{fo). 

Substituting ([530]) and (f53T]) into ([S^lQ]) . we obtain 



pI(x) > e-^+^^ inf y f^{x, - C2{v)Tf, v)dv 

> e-^+^^ inf ( / fm{xi - vTf,v)dv - CNg{h){2 + Tf){Ax + Av) 

1 ^IEl. f 

> e '^+^tmf foR{xi-vTf,v)dv. 

In the last line, we used the smallness assumption (|5.45|) to see 

iVg(/o)(2 + Tf)iAx + Av)<^f Mix - vTf,v)dv. 
Hence we have from ()3.27p 

We now turn to the proof of the lower bound of T^. Note from the estimate (|5.32p of 
Lemma [521 we have 



PR 



Now the pointwise upper bound estimate of U and T follows directly from the following 
observations: 

^n. ^ < Jfj^{x,v)\v-U^{xrdv 

n 

The following continuity property of local Maxwellians is from |15j . 

Lemma 5.6. |15j Let f and g be solutions of il.l\) such that Np{f) < oo and Np{g) < oo. 
Suppose f and g satisfy estimates ()3.28p . Then we have 

(5.52) \\M{f) - M{g)\\L^^<CT,\\f - g\\Li- 

The proof of this lemma is given in [15], we present here the detailed proof for the reader's 
convenience. 
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Proof. We observe from Taylor's theorem 

M{f)-M{g) = M{pf,Uf,Tf)-M{pg,Ug,Tg) 

dM{e) 



dM{e) ,, dM{e) 



where 



dp 

dMje) _ dM 
dX " 'dX 



dUi 



T 



f-g- 



dT 



9/+(l-e)g 



Multiplymg (1 + \v\'^) to both sides and mtegrating with respect to {x,v), we get 



'XT'* 



\Mif)-Mig)\il + \v\'^) dvdx 



(5.53) 



Pf^g 

+ I Uf.g 

+ I Tf^g 



dM{e) 



dp 
dM{e) 



dUi 
dM{9) 



dT 



(1 + \v\'^)dvdx 

(1 + \v\ )dvdx 
(1 + |t;| )dvdx. 



We then substitute the following estimates into (|5.53p 



dM 

dp 

dM 

dM 

dT 



1 



\v-ur 



-e 2T 



V(27rT)^ ^ T 
d 1 p 



)e 2T , 



to obtain 



/, 



dM{e) 



dp 



2ry(2^ 



{l + \v\'^)dv 



-e 2T -\- 



e 2T 



7(2^^ 2r2 



Pe 



\^-Ue\ 



2T, 



Pe jRrf yJ{2TTTg 



N 



e {l + \v\'')dv 



Pe 



l^-upr 



pe jRd y/{2TrTe)^ 



e 2^e dv 



+ - 



Pe 



-\^-Uer 



pe jRd y/{27rTe)^ 



\v\ dv 



= —pe + —ipe\Ue\'^ + dpeTe) 
Pe Pe 

= l + \Ue\^ + dTg. 
In a similar manner, we can estimate the remaining terms of (j5.53p as follows: 



dMie) 



dU 
dM{d) 



{l + \v\'^)dv < 



dT 



^e 
Pe . 



[l + \Ue\'' + Te), 



{l + \v\')dv < ^{l + \Ue\' + Te) 
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Substituting these estimates into (|5.53p . we get 

\Mif)-M{gm + \v\^)dv 
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(5.54) 



= {l + \Ue\^ + Tg)pf_g + ^(1 + \Ue\' + Te)\Uf.g\ 
J-ft 



We then integrate over T^; and employ Lemma 15.21 to have 

(5.55) \\MU) - M{g)\\Li <C f (pj^g + \Uf.g\ + Tf.g)dx, 

where C = C(iVg(/o), e"*). We estimate these terms separately. Note first that we have 

(5.56) / pf-gdx< \f -g\dvdx <\\f - gW^i. 
On the other hand, we observe that 

\Uf-g\dx = C{t) [ \Uf-Ug\\pf\dx {■.■ p>Ce~^'^) 
Jt3 



Ct / \UfPf - UgPf\dx 



Ct / \UfPf - UgPg + UgPg " UgPf\dX 

Ct \UfPf-UgPg\dx + CT \Ug\\pg - pf\dx. 



The first term can be estimated as follows 



/ \UfPf -UgPg\dx = fvdv- / gvdv 

JTd Jfd J^i J^d 

1/ — g\\v\dvdx 



dx 



< 

< Wf-ghi- 



Similarly, we have 



\Ug\\pg - pf\dx < C{t) / 1/ - g\{l + \v\')dvdx. 

Combining these estimates, we get 

(5.57) f \Uf^g\dx<CT\\f-g\\Ll. 

Finally, we consider 



'S-9 



dx 



C{t) I \Tf-Tg\\pf\dx (•.• p>Ce-*) 
/ \TfPf-TgPf\dx 

/ \TfPf - TgPg + TgPg - TgPf\dX 
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/ \TfPf-TgPg\dx+ \Tg\\pg- pf\dx. 



The first term can be estimated as follows 

dx 



\TfPf -'^gPgldx =0/ / f\v?dv-l g\v\'^dv 

< - / \f -g\\v\dvdx 

< cf \f-g\{l + \v\'^)dvdx. 
Similarly, we have 



\Tg\\pg - pf\dx < C{t) / 1/ - g\{l + \v\')dvdx. 



The combination of these two estimates yields 

(5.58) / \Tf_g\dx<CT\\f-g\\Ll- 

We now substitute dESSD, dSSZD, (fH35D into (153^ to obtain 

\\M{f)-M{g)\\n<CT\\h-g4L.. 



U 



The above continuity property also holds for discrete distribution functions. 

Lemma 5.7. Let f^ and g^ numerical solutions defined recursively bv \2.26\ corresponding 
to initial discretization /^ and g^ respectively. Assume that iVq(/^) < cxo and Nq{g^) < 00 
with g > 4. Then we have 

(5.59) \\M{f^)-M{gl)h^<CT,m-gMLl■ 

Proof. Combination of Lemma 15.41 Lemma [5.51 and Lemma lS. 61 gives the desired result. D 

6. Consistent form 

In this section, we transform the BGK equation (jl.ip to derive a consistent form which 
is compatible with the reformulated scheme (j2.26p . To be consistent with the notation for 
discrete problems, we employ the following notation for the smooth distrbution function: 

f{x,v,t) = f{x-vAt,v,t). 

Theorem 6.1. Let f be a smooth solution of (jl.ip . then we have 



(6.60) 
where 



f{x,v,t + At) = ^—f{x,v,t) + ^—M{f){x,v,s) 

K + At K + At 

At _ 1 _ 

+ 7-^1 + ^■^2- 

K + At K + At 



7^l = MU)ix,v,s)-M{f){x,v,s) 
T^2 = Rm — Rfi 
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and 



r-t+At ^J^ 

Rm = {s-t){-j^+vVM){xe^,v,te^) 

/t+At 
(s - t - At)Mif){xe2,v,te2)ds. 

Proof. Along the characteristic hne, we have from (ll.ip 

df 1 

-^(x + vt, V, t) = -{Mf - f){x + vt, V, t). 
at K 

We integrate in time from t to t + At to obtain 

I rt+At 

f{x + {t + At)v, V, t + At) = f{x + tv, v,t) + - / {Mf - f){x + vs, v, s)ds, 

1^ Jt 

or, equivalently, 

(6.61) f{x,v,t + At) 

I rt+At 

= f{x - Atv, v,t) + - ( Mf - f ){x - {t + At - s)v, V, s)ds 

1^ Jt 

= f{x-Atv,v,t) + -{lM+If)- 

K 

AppUcation of Taylor's theorem around {x — Atv, v, t) gives 

M{x- {t + At- s)v,v,s) = M{x- Atv,v,t) 

+ {s-t)vV^{M{x0^,v,te^)) 

+ {S-t)—^{XQ^,V,tQ^). 

where xq^ lies between x — Atv and x — {t + At — s) and tg-^ £ [s, t]. Hence we have 

(6.62) Im = AtM{x- Atv,v,t) + RM■ 
wh.ere 

/t+At 
{s - t)v ■VxiM{xg^,V,td^)) 

f^+^\ ,dM, 

+ / {s-t)—^[xB^,V,t0^). 

On the other hand, we employ Taylor's theorem around (x, v,t + At) to get 

f{x-{t + At-s)v,v,s) = f{x,v,t + At) 

+ {s-t- At)v-V;,{f{xe2,v,tg^)) 

. .df 
+ [s-t-At) — {xe2,v,tg2), 

where xg^ lies between x and x — (t + At — s) and tg^ E [s, t]. This gives 

(6.63) If = Atf{x,v,t + At) + Rf, 
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where 

/t+At 
{s-t-At)vVifixe„v,te,)) 

rt+At ^f 



+ / {s-t-At) — {x02,v,t02) 

r-t+At 

{s - t - /\t)M{f){xe.^,v,te^). 



n 
for some appropriate 62 G [t-,t + At]. Substituting (|6.62p and (|6.63p into (j6.6ip . we obtain 

f{x,v,t + /\t) = f{x- Atv,v,t)^ M{x- Atv,v,s) f{x,v,t + At) 

+ -{Rm-Rj). 

K 

We then cohect relevant terms to have 

/(i, t), ( + A() = ;rj^/(^' - Ato, !>, () + ;-^>l (f){x - Ato, o, s) 

f{x, V, t) + ■M(/)(x, V, s) 



K + At ^ ' ' ' K + At 
.64) + ^( i?^ - R,) 



K ~ At ~ 

-f{x, V, t) + M (/)(x, V, s) 



K + Ar ' ' K + At 
+ ^^^^[Mif){x,v,s) - M{f)ix,v,s)) 

{Rm- Rf)- 



K + At 

Finally, we put 

7^l = M{f){x,v,s)-M{f){x,v,s), 
^2 = Rm ~ Rf- 

to obtain the desired result. D 

Before we estimate these remainder terms, we need to establish the following technical 
lemmas. 

Lemma 6.1. Let f be a smooth solution of l\l.l\) corresponding to an initial data /q. Then 
we have for q > d+ 2 

\\j^ip,U,T)\\Li + \\V,ip,U,T)\\Li + \\v ■V,{p,U,T)hi <C,N,{fo). 
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Proof. We prove only the first estimate. Other estimates can be treated in a similar manner. 

For ^{x,v) = l,v, \v — U{x)\'^, we have 



(6.65) 
and 



f{x,v,tMv)dv <iV,(/) 



— / f{x,v,t)ip{v)dv 



"'"I ■ 'l + \v\)i 
< C.Ngifo) 



— {x,v,t)ip{v)dv 



dv 



{v-Vf+'^{M-f))v{v)dv 



.66) 



Therefore, we have 
(6.67) 



V ■ V fip{v)dv 



- '^'Vr^' '(1 + 1^1)'^ 



p\, \pU\, \pT\ < C,N,{fo), 
3p OipU) OipT) ^ ^^^^^^^^_ 



dt 



dt 



dt 



Then the result follows from (j6.67p and the lower bound estimates for local density and 
temperature given in Theorem 13.11 D 



Lemma 6.2. Let f be a smooth solution corresponding to an initial data /q. Then we have 
for q>d + 2 

II — 



|m + llV^A^II^i + \\vV^,M\\li ^ C.^iifo)- 



Proof. We recall the chain rule: 
dM 



dpdM dU dM dT dM 

If ~ 'dt~di^~dt"dlA^~dt~dT' 

^ ,, ^ dM ^ ^^ dM ^ dM 

VxM = V^p^— + V^U-—- + V^T—-, 
dp dUi dT 

^ . . ^ dM ^ ^dM ^ dM 

V^M = v-V^p—- + v-V^U--— + vV^T—- 

dp dUi dT 

and apply the estimates of the Lemma 16.11 to obtain 

„dM 



dt 



\l}, + WxM\\li + \\v-VxM\\li 



< CgN.ifo) 



dM 



dp 



+ 



dM 



dUi 



+ 



dM 



dT 



{l + \v\fdv 



< C,N,{fo) [{1 + |[/|2 + T) + -^(1 + |C/|2 + D 

< C,N,{fo). 



This completes the proof. 



D 
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The following lemma provides the estimate of the remainder term, which will be crucial 
for the convergence proof. 

Proposition 6.1. i?x, Rf satisfies following estimates. 

(6.68) ||7^l||il < CT,||/-7||^i + C,iV,(/o)Ai, 

(6.69) ||7^2||il < C,N,ifo)At^ 

Proof. We observe that 

11/ - /ILi < \\Atv . VJWli < C,Ngif)At. 
This, with the estimate of Lemma |6.2|, yields 



^2 



l^ilLi = ||A^(/)-A^(/)ILi 

< \\M{f)-M{f)\\Li + \\M{f) - M{f)\\Li 

< CT,\\f-f\\Li + \\AtvV,M{f)\\L^ 

< Ct, 11/ -7lLi+C,iV,(/o)Ai. 



On the other hand, we have by Lemma 167 

rt+At 



ll^2|L^ < J^ {s-t)[\\vVM{f)hl + \\j^M{f)h^ + \\M{f)h.)ds 
< CgNg{fo)l {s-t)ds 



< C,iV,(/o)(At)2. 

D 

7. Proof of the main result. 

Before we delve into the proof of the main theorem, we need to establish some technical 
lemmas. 

Lemma 7.1. Suppose Ng{fQ), Ng{f^) < oo with q > d + 2, then we have 
\\f- IrWlI ^ 11^ - fnU + Cg{Ng{fo) + iV,(/0)) (Ax + AvAt). 
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Proof. We divide the estimate of ||/ — //jH^i into the following four integrals. 

\\f - fliWLi, 

= I \f{x - vAt,v,t) - f^{x,v)\{l + \v\fdxdv 

< / \f{x-vAt,v,t) - f{x - C2iv)At,v,t)\{l + \v\fdxdv 

(7-70) + / |/(x - C2{v)At,v,t) - fl{x - C2{v)At,v)\{l + \v\fdxdv 

+ / \fl{x - C2{v)At,v,t) - mCiix) - C2{v)At,v)\{l + \v\)^dxdv 

+ / IfniCiix) - C2iv)At,v,t) - mx,v)\{l + \v\fdxdv 
= 1 + 11 + 111 + IV. 
By Taylor's theorem, / can be treated as follows 

/ = I \f{x-vAt,v,t)- f{x-C2{v)At,v,t)\{l + \v\fdxdv 
\v - C2{v)\At\d^f{xe,v,t)\{l + \v\fdxdv 

< AvAt \d^f{xg,v,t)\{l + \v\fdxdv {■ .• \v - C2{v)\ < Av) 

< AvAtN.if) J j^^^{l + \v\fdxdv 

= CAvAtNgifo). 
On the other hand, we have by the change of variables with respect to x 

II = I \f{x-C2{v),v,t)-fl{x-C2{v)At,v)\{l + \v\fdxdv 

\f{x,v,t) - fl{x,v)\{l + \v\fdxdv 

= \\f - IrWlI- 
We now turn to the estimate of ///. We define for simplicity Rij as 

/// = I \fl{x-C2{v)At,v)-fl{Ci{x)-C2{v)At,v)\{l + \v\fdxdv 
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/I 



Xi 



x — VjAt — Xs 
Ax 



'J's,j,R + 



Xs — {x — VjAt) 
A^ 



s E 



VjAt-Xs Xs+i - {xj -VjAt) 

~Ax ^^+i'^'^ + 

X-VjAt-Xs^ ri 
Js.i 



Ax ^^■^■'^. 



rn 
Js-l,j,Rj 

(1 + \v\)'^A^j{x,v)dxdv 



Ax 



Xi — VjAt — x. 



S fU 



Ax 



fs+l,j,R + 



Xs-ix-VjAt) 

•^■'^^ A^ ^^-i'^'^. 

Xs+i - {xi -VjAt) 



Ax 



Js,j, 



R 



(1 + \v\)'^Af^j{x,v)dxdv 
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- H/ Rij {'^ + \'"\f-^i^j{x,v)dxdv. 

We first observe that if 2; G [xi_i,Xj) and Xi — VjAt G [xs,Xs+i), tfien we have either 
X — VjAt £ [xs-i,Xs) ov X — VjAt E [xs,Xs+i). Hence we divide the estimate into the 
following two cases, 
(i) The case of x — VjAt £ [xs-i,Xs): For brevity, we put 

X — VjAt — Xs-i 
Ax 



5. . = Zl ^^_ 

«J — A 

Ax 



to see 



^d — I (^i,jfs,j,R + (1 ~ 0-i,j)fs-l,j,R - ( Kjfs+l,j,R + (1 ~ h,j)fsJ,R ) 
— ^i,j\fs+l,j,R ~ /"j,ijl + (1 ~ 0,i,j)\fs,jR ~ fs-l,j,R\ 

(ii) The case of x — VjAt G [xg, Xg^i): An almost identical argument gives 

Ax 



(1 + l^ii 



Rii < 2N„(fV,) — 
From (i) and (ii), we have 

/ Rij {l + \v\fA^j{x,v)dxdv 



^^'^'(^«'r?r'(T^"^i"i>^''^''" 






- ''^''^«'/"? C (T^" + 1"' + 5^""''""" 



VI f^^+l Ax 



which gives 



^.r(fn. C,{AxrAv 



< C,iV,(/|^)Ax. 



Finally we estimate IV. Suppose Xj — VjAt G [xs,Xs+i), which implies Xj+i — WjAt S 
[xs+i,Xs+2)- This gives for {x,v) G [xj_i,Xj) x [v-_i,v-,i) 

J 2 J^ 2 

fu,, „.A. „^_^i^^^j^i^^.n I x.+i-(x,-z;,At) 

jR\Xi - VjlXt,V) - — Js+l,j,R H -^ JsJ,R 



CONVERGENCE OF A SEMI-LAGRANGIAN SCHEME 



29 



and 



m^^v) 



Ax ^^+1'^'^+ Ax ^^'^'^^ 

X- Xi fXi+l - VjAt - Xs+l „ 



+ 



Ax V Ax 

X^-j-l — X / Xi — VjL-\t — Xc 



fs 
/s+1, 



Xs+2 - {xi+1 - VjAt) 



+-2,j,R 



Ax V Ax •'-■+^.i.« 

Hence we have for (x,u) E [xj_i,x.j) x \v-_i^v-,i) 

■I ■?. J'^ -7. 



+ 



Ax 
Xg+i - {xj - VjAt) 

Ax 



fl 



+1J,R 



rn 



f^iCiix)-C2{v),v)-f^{x,v) 

X- Xi (Xi+i-VjAt-Xs+l n , Xs+2.R - {xi+l - Vj^t) ^„ 

Js+2,j,R H TT- Js+l,j,R 



Ax 



Ax 



'^^ V a;; ■'^+ij,i? 



' J s 



Ax 
, x^+i - [xi - VjAt) ^„ 



Therefore, we have by a similar argument as in /// 



IV < C,N,Ul)Y^ 






Ax 



Ax 



(Axl^Au 



»j 



{-^ + \vo\Y-^ 



< C,N,{fl)Ax, 



where we used 
result. 



Ax 



< 1. Substituting these estimates into (j7.70p . we obtain the desired 

D 



Lemma 7.2. Suppose Nq{fQji), Nq{ff^) < oo with q > d + 2, then we have 

(a) \\M{Tr) -Mim^i < Cr.Wf - m^ + C,N,{foR){Ax + Az;At), 

(b) \\Mif^)-M-if^)hi < C,N,ifo)Av, 

(c) \\M-if^)-E{M^f^))h^<C,N,ifo)iAx + Av). 

Proof. By Lemma 15.61 and Lemma 17.11 (o) , we have 

WMif) - M{f^)h.^ < CTf\\f-mq 



< Ct, 



mLl + C,N,{fo){Ax + AvAt), 



which concludes (a). 
We now prove (b). Let 

Then we have 

M{f^)-M-{f^) 

f^R 



Vi = {^^i_l < \v\ < V. i}. 



(27rrg)^ 



■ exp(- 



- y^i'^ -c'2(^^)isup 

. Vi 



111 _ r/'"!^ 
P ^r\ ■ 

2T^ ' 

dMiJl) 



E 



Pr 



■ exp( 



(2vrT]J) 



"n\N 



2T\ 



)Xv. 



R 



dv 



Xvi 
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< 
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Av y^ sup 



dv 



Xvi- 



To estimate this, we recall Lemma 15.51 to see 









dv 


< 


PR be 

^{2TrT^Y V 


'-'rI 

' ^R 


exp 


( 


\ve 


-ui 


|2 






2TI 












< 


n ( 1^^ 
Ct exp 1 


-VI 
2TI 


2x 


















< 


CTexp(— Ct^) 


w 












where 


we 


used 


























( be 

exp 

V 


2TI 


-) = exp ( 


\Vi - 


- {Vt 


-Vq) - 

2TI 


-^l? 


■) 



< 



exp 



2\{vi-ve)^Vl 



n\2 



2T^ 



Vi 



^ ^"Pl-^J^"Pl ^ 



R 



\Av\^ + \U^\'^ 



R 



Hence we have 



\\M{f^) - -M"(7|^)|Li < C,N,ifo)Av. 



We now turn to (c), We consider 

\\M-m-E{M^{f^))\\^^ 



Ell^^(-^"(/^)(^'^)--^"(/^)(^-+i'^^))<^ 



ml 



Ax 



+E 



Xm+l ^ 



m,i 

A + B. 



Ax 



M-{mx,v)-M-{mx^,ve)]Ai^^, 



Lh 



LI 



Since ^^-^ < 1, we have 



L\ 



A < Y.\\M'^{f^Xx,v)-M-{f^){x,„.^i,vMm,i 

= E/, \M''{mx,v)-M''{f^){xm+i,ve)\{l + \v\fdvdx 

m,£ -^^m.i 

/A^(x v) 



+ \v\)i-^ 



dxdv 



< CgNgifo)iAx + Av). 
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Here we used for {x,v) £ [xm,Xm+i) x [v^_i,v^^i) 

M^imx,v)-M^{mxm+uv, 

< 
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(x - x^){d^p + d^U + d^T) ■ (^— + ^j77 + -T^) + {v- Vi)- 



< 



N,iW 



dp 



dU 



dT 



dv 



(A3; + At>), 



il + \v\)i-^ 

which can be obtained by an almost identical argument as the one given in the proof of 
Lemma 15.31 B can be handled in a similar manner. 



D 



We are now in a position to prove our main theorem. 



7.1. Proof of the main theorem. We subtract (j2.26p from (j6.64p and take L2 norms to 
obtain 

\Ll = —^M-,nAt)-fU;-)\\Ll 



||/(.,.,(n + l)At)-/^+i(v 



(7.71) 



K + At 

K + At" "^2 K + At" "^2' 



where 

TZi = M{f){x,v,s)-M{f){x,v,s), 
^2 = Rm ~ Rf- 
We then collect the estimates of Lemma 17.11 and Lemma 17.21 

WI^rWlI < \\f - mLi+C,N,{fo)AvAt^+CgN,{fo)Ax, 
\\M{f)-E{M{fm,l < \\M{fl-M{n)h. + \\M{f^)-M^{f^)h. 

+ WM^if^) - E({M-{f^)))h^ 

< C\\f- fjhi + C,Ng{fo)AvAt + C,iV,(/o)(A^ + Ax) 
||7^l||il < C\\f-f\\Li+C,N,{fo)At, 
\n2h1 ^ CN,{fo){Atf. 

We substitute these estimates into (|7.7ip to obtain 

||/(-,-,(n + l)At)-/^+i||M 



< 



11/ - fRhl + C,,TfN,{fo)AvAt + C,,TiVg(/o)Ax 



K + At 

+TT^( 11/ - IrWlI + N,{fo)AvAt + N,{fo){Av + Ax) 

Cg,T^At 



K + At 



( 11/ - IrIIli + ^^(/o)^0 + ^''^«(^o)^Ta^(^*)' 



CT,At\ 
^ + ^h^) f-fn 

K + At/ 

1+ \ )At;At 

K + At 
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CT,At\ 

1 + ;^^)^^^ + ^^) 



^''^^ {Atf. 



K + At 
We now put T = jf . , for simplicity of notation and iterate the above inequality to obtain 



K+At 



||/(.,.,iv,At)-/^'||^. <(i + r)^*||/(.,.,o)-/° 11^1 

Nt 



(7.72) ^* 



+ ^(l + r)^(Ax + Ai 



Aft-1 



Note that from the elementary inequality (1 + x)" < e"^, we have 

(1 + r)^* < e^*r 

< Q—mrET- 

= e ''+^* , 



and 



i=l 



(i + r)-i 

^ ^ f£L>K + At 

< CT,e «+^* — . 

f At 



Similarly, we have 



5:(i + r)'<c.,e^^. 

i=0 



We then substitute these estimates into (|7.72p to obtain 



\\f{; ■,NtAt) - f^'hi < e^^W II /(., .,0) - /: 







-2 - ■■ ■'-R 11^2 



+ CTfe -+A* (k + At) Ax 

^^•^^^ ^ ^^>^ , ,,(Ax + At;) 

+ CT^e-+A* {K + Aty ' 



At 



+ CTfe -+A* At. 



To estimate ||/o — /{^IIlI) we decompose the integral as 

ll/o - //?IIli ^ Wfo- fR\\^(H<R) + \\fo\\m\v\>R)- 
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By direct estimates, we have 

Wfo- fR\\Ll(\v\<R) = ll/o - -E'(/0'^|i;|<_r)IIl1 

< C,Ng{fo){Ax + Av). 
On the other hand, the second term can be estimated as fohows 

II/oIIli{H>r) < Ngifo) I _ „,... dv 



|>ij(l + bP'' 



< C,N,{h) 



We substitute them into (j7.73p and take 

C{Tf,q,fo)= max{CT,e^,C7,iV,(/o)}, 
to obtain 

\\f{.,.,Tf)-f^^\\^.<C{Tj,qJo)(Ax + Av + ^^^±^ + At + 



^2 



At {1 + R)i+^ 



This completes the proof. 
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